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i—i The work reported below is a first of its kind study of the properties of turbulent flow without 

^ strong mean shear in a Newtonian fluid in proximity of the turbulent /non-turbulent interface, 

o 

with emphasis on the small scale aspects. The main tools used are a three-dimensional particle 

in 

o 

tracking system (3D-PTV) allowing to measure and follow in a Lagrangian manner the field of 



X 



velocity derivatives and direct numerical simulations (DNS). The comparison of flow properties 
in the turbulent (A), intermediate (B) and non-turbulent (C) regions in the proximity of the 
interface allows for direct observation of the key physical processes underlying the entrainment 
phenomenon. The differences between small scale strain and enstrophy are striking and point to 
the definite scenario of turbulent entrainment via the viscous forces originating in strain. 
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Turbulent entrainment (TE) is a process of continuous transitions from laminar to tur- 
bulent flow through the boundary (hereafter referred to as interface) between the two co- 
existing regions of laminar and turbulent state. This process is one of the most ubiquitous 
phenomena in nature and technology since, in fact, most turbulent flows are partly turbu- 
lent: boundary layers, all free shear turbulent flows (jets, plumes, wakes, mixing layers), 
penetrative convection and mixing layers in the atmosphere and ocean, gravity currents, 
avalanches and clear-air turbulence. 

The first physically qualitative distinction between turbulent and non-turbulent regions, 
made by CorrsirP^, is that turbulent regions are rotational, whereas the non-turbulent ones 
are (practically) potential, thus employing one of the main differences between turbulent 
flow and its random irrotational counterpart on the 'other' side of the interface separating 
them. 

The main mechanism by which non-turbulent fluid becomes turbulent as it 'crosses' the 
interface is believed to involve viscous diffusion of vorticity {yujiV 2 uJi) at the interface^. 
Corrsin and Kistler 2 also conjectured that the stretching of vortex lines in the presence 
of a local gradient in vorticity at the interface leads to a steepening of this gradient since 
the rate of production of vorticity is proportional to the vorticity present. The mentioned 
processes are associated with the small scales of the flow. However, at large Reynolds 
numbers the entrainment rate and the propagation velocity of the interface relative to the 
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fluid are known to be independent of viscosity (see Ref. [3B1I5] for more information and 
references). Therefore the slow process of diffusion into the ambient fluid must be accelerated 
by interaction with velocity fields of eddies of all sizes, from viscous eddies to the energy- 
containing eddies, so that the overall rate of entrainment is set by the large-scale parameters 
of the flow. This means that although the spreading is brought about by small eddies 
(viscosity), its rate is governed by larger eddies. The total area of the interface, over which 
the spreading is occurring at any instant, is determined by these larger eddies^. 

Until recently it was difficult to implement Corrsin's distinction: it requires information 
on small scale vorticity and strain which experimentally was not accessible. This is why very 
little is known about the processes at small scales and in the proximity of the entrainment 
interface. A few exceptions are recent particle image velocimetry (PIV) and planar laser- 
induced fluorescence (PLIF) experiments by Westerweel et alP^of a jet and experiments by 
Holzner et al. 8 on oscillating grid turbulence, in addition to the direct numerical simulations 
(DNS) of a temporally developing plane wake by Bisset et alP and a similar analysis of 
an axisymmetric configuration by Mathew and Basu 1Q . Unfortunately, in the PIV/PLIF 
experiments only the azimuthal component of vorticity in a two-dimensional cross-section 
was accessible. 

The main objective of the study presented is a systematic analysis of the small scale 
dynamics associated with turbulent entrainment. The special emphasis is on the processes 
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involving the field of vorticity, u>i, and its production/destruction by inertial UiUjSij and 
viscous vuJiW 2 Ui processes in the proximity of the interface (s^ are the components of the 
fluctuating rate-of-strain tensor, v is the kinematic viscosity). Studying the production of 
vorticity requires access to the field of strain as well (and thereby also to the dissipation, 
2vs.ijSij). We briefly recall the local budget equations for enstrophy, reading = tOiUjjSij + 
uUiV 2 Ui and strain, written as -^y = -SijSjkSki - \uiUjSij - a fj x . + //.s (J V 2 .s, r 

Experimentally, a turbulent / non-turbulent interface was realized by using the oscillating 
planar grid described in Holzner et alP. The grid is a fine woven screen installed near the 
upper edge of a water filled glass tank, oscillating at a frequency of 6 Hz and an amplitude of 
4 mm. The scanning method of 3D particle tracking velocimetry (3D-PTV) used here is pre- 
sented in Hoyer et al.^. In order to access the viscous term of vorticity, the postprocessing 
procedure was extended. Due to experimental noise it is not possible to obtain the Lapla- 
cian of vorticity, z/V 2 o;, directly through differentiation as it involves second derivatives of 
the vorticity field. Instead, the viscous term is obtained from the local balance equation of 
vorticity in the form V x a = uV 2 u: by evaluating the term V x a from the Lagrangian 
tracking data. The derivatives of Lagrangian acceleration, ddi/dxj, are calculated in the 
same way as the derivatives of the velocity, dui/dxj, described in Hoyer et al.^ and subse- 
quently they are interpolated on a Eulerian grid. The number of tracked particles is about 
6-10 3 in a volume of 2x2x1.5 cm 3 and the interparticle distance is about 1 mm, which is 
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slightly above the estimated Kolmogorov length scale, 77=0. 6mm. The Taylor microscale, 
A, is about 7 mm. The spacing of the Eulerian grid was taken equal to the interparticle 
distance. Further details on the data analysis will be discussed in the forthcoming full paper. 
In both experiment and simulation, the Taylor microscale Reynolds number is ReA=50. 

Direct numerical simulation (DNS) was performed in a box (side-lengths L ± , L 2 , L 3 ) 
of initially still fluid. Random (in space and time) velocity perturbations are applied at 
the boundary x 2 =0. The procedure of generating the boundary conditions is as follows. 
For a fixed time and in the discrete set of points, X\ = kA h x 3 = lA t (k, 1 - integers), 
each velocity component, Ui (i = 1, 2, 3), is calculated as Ui = V^, where £ is a random 
number within the interval [—1,1] and Vi is a given velocity amplitude. For other times and 
spatial points (xi, #3) boundary velocities are obtained by cubic interpolation in time and 
bilinear interpolation in space. At each time the three boundary velocity components yield 
zero average value over the boundary plane. The method of boundary velocity assignment 
determines the velocity scale, V = max(Vi) and the length scale A;. Together with the 
viscosity of a fluid, v, these parameters define the Reynolds number Re = VAi/u = 1000 of 
the simulation. The Navier-Stokes equations are solved with periodic boundary conditions 
for the directions x 1 and x 3 , with periods L 1 and L 3 , respectively. The computational domain 
is finite in the x 2 direction, as x 2 < L 2 . Shear-free conditions du\jdx 2 = du 3 /dx 2 = u 2 = 
are imposed at the boundary x 2 = L 2 . A mixed spectral-finite-difference method is used for 
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the spatial discretization and the time advancement is computed by a semi-implicit Runge- 
Kutta method^. The resolution is 192x192 Fourier modes in x± and 23 directions and 192 
grid points in X2 direction. The local Kolmogorov length scale is twice the grid spacing. 




FIG. 1: Average profiles of uj 2 {x2) and 2s 2 (x2) from PTV (symbols), and DNS (lines) relative to the 
interface, X2 = X2 — x^ on linear scale (left) and on log-scale (right). The values are normalized by 
the respective maxima (max(w 2 )~ max(2s 2 ) ~3.6 s -2 for PTV and 3.0 for DNS). The horizontal 
error bars display the sensitivity of the experimental result on the enstrophy threshold, 0.5-2.5 
s~ 2 , the vertical error bars represent the accuracy of the measurement. The symbol (•) denotes 
the ensemble average of the respective quantity 

In both experiment and simulation, turbulence is generated at the plane x<i =0 and 

propagates along x 2 > 0. Firstly, the interface is identified at x\{t\ using a fixed threshold 

of enstrophy, (for details see Ref. [8] and references therein) and the analysis is done with 

respect to the interface location, as in Fig. 1, in which profiles of enstrophy, u 2 = UiUi 

and strain rate, s 2 = s^s^-, averaged over homogeneous x±,x 3 directions on linear scale 

(left) and log scale (right), are shown. The distance to the interface, 22 is normalized by 
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the Kolmogorov length scale, r\. The 'proximity' or 'region of the interface' hereafter refers 
to the interval -5< £2/77 <5. We observe that the rate of strain on the non-turbulent 
side of the interface remains high in contrast to enstrophy which drops much more steeply. 
Experimentally it is not possible to obtain enstrophy lower than a small (but finite) level of 
noise. This is one of the reasons why comparison to DNS is presented for all the results. In 
the DNS, the numerical noise level is reached at £2/?] >10 and this level is about 25 decades 
lower in magnitude, see the inset in Fig. 1 (right). 




£ 2 h x 2 h 



FIG. 2: Average profiles of production and viscous destruction terms of strain and enstrophy (left) 
and their rates (right). The inset shows the individual Lagrangian trajectories of z^cujV 2 ^ obtained 
from PTV. Lines are from DNS, symbols are from PTV 

Fig. 2 shows profiles of production and viscous terms of strain and enstrophy (left) and 
their rates (right). We note that the viscous term, z/tUjV 2 ^, exhibits a remarkable behavior 
showing a distinct maximum in the region of the interface. In addition, the individual 
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Lagrangian trajectories (examples are shown in the inset in Fig. 2, where the abscissa is 
i — t — f , normalized by the Kolmogorov time scale, t v , and t* is analogous to x* 2 ) possess 
such an extremum. Therefore, we use the maximum of the viscous term as the exact location 
of the interface, defined in a physically more appealing way than the threshold-dependent 
crossing of x 2 /r]=0. For the further analysis we define three physically distinct regions of 
the interface with respect to the maximum of vu),jS/ 2 uJi (marked in Fig. 2): (A) the turbulent 
region, in which the behavior of the viscous term is 'normal', i.e. it is negative in the mean, 
(B) the interval between the peak and the point where (z/a;jV 2 cjj)=0 is termed intermediate 
region (with the 'abnormal' viscous production) and, (C) the non-turbulent region from the 
peak to rr 2 /?7=5. The positiveness of both UiUjSij and vu)iV 2 u)i is a peculiar feature of the 
regions B and C, in contrast to region A, where, in the mean, z/WjV 2 u;j contributes to the 
destruction and UiUjSij to the production of uj 2 . It is noteworthy that strain behaves rather 
differently from vorticity. In particular, the viscous term, z/s^-V 2 ,^, is negative in the mean 
in all three regions, i.e. it is not building up s 2 . In Fig. 2 we see also that strain production, 
-SijSjkSki, is significant and it is (in the mean) not balanced by usijV 2 Sij in region C. When 
the rates of quantities are considered it appears (Fig. 2, right) that the role of viscous 
production is even more important: the term voj^uJijiJ 1 attains high positive values in 
region C, decreases along region B and finally becomes negative in region A (balancing the 
average UiUjSij/u 2 ). In contrast, the term UiUjSij/u 2 and the analogous rates of the strain 
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viscous and production terms do not change as drastically and remain of the same sign. 
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FIG. 3: PDF of various quantities from experiment (top) and simulation (bottom) according to 
the division in 3 regions: turbulent (left), intermediate (center) and non-turbulent (right). UiUJjSij 
(-), WjV 2 w, (- -), -SijSjkSki (- ■), vsijV 2 Sij (• • • , only bottom) 



Fig. 3 presents the estimations of probability density functions (PDFs) of the relevant 
terms from the different regions (A,B, and C, from left to right; PTV top, DNS bottom). 
Consistently with the other results, the PDFs of both UiUjSij and vuji¥ 2 uji are positively 
skewed in regions B and C. In region B we note that the probability of negative events 
of vuJiW 2 uJi and positive events of UiUjSij increases as compared to region C. Finally, as 
expected, in region A the PDF of z/tUjV 2 ^ is negatively skewed. The changes of the strain 
production and viscous terms between the regions A-C are less drastic. Essentially, -SijSjkSki 
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is positively and z/SjjV 2 Sjj is negatively skewed in all the three regions. 
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FIG. 4: (left) Cosines of the angle between vorticity and W (-,o) and between vorticity and vX7 2 u> 
(- -,*). (right) Invariants, Q (-, o), R (- -,*), S (— • — , only DNS). Lines are from DNS, symbols 
are from PTV 

For the understanding of the interaction of strain and enstrophy in the proximity of the 



interface it is very instructive to look at the invariants of the gradient tensor: Q = lo 2 — 2s 2 



— V 2 P shows the relative strength of vorticity and strain, R=-l/3(sijSjkSki +3/Au>iU)jSij ) 
relates to their production terms, and S is the quantity related to the two viscous terms, 
S = vujiV 2 0Ji — 2vs i j'V 2 Sij (Fig. 4, right). Since the mean values of Q,R and S vanish 
identically for homogeneous turbulence, their nonzero values indicate the degree of inhomo- 
geneity in the proximity of the interface. Apparently, inhomogeneity is the property which 
is maximal where also vuJiV 2 uJi is maximal. In the same context it is also interesting to look 
at the cosine of the angle between vorticity and its Laplacian, V 2 ^, shown in Fig. 4 (left), 
which exhibits significant changes across the regions A,B, and C. The observed transition 
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from positive (alignment) to negative (anti-alignment) values is in agreement with the qual- 
itatively different behavior of z/u;iV 2 u;i in these regions. In contrast to that, the alignment 
between vorticity and the vortex stretching vector, Wi = ujjSij, changes only weakly through- 
out regions A-C and is consistent with the positiveness of (uJiUjSij) mentioned above. The 
results indicate that an interpretation of the viscous term i/cUjV 2 ^ as interaction between 
strain and vorticity due to viscosity (i.e. due to the curl of the viscous force originating from 
the divergence of the strain tensor) is physically more appealing than 'simple' diffusion of 
vorticity due to viscosity. We emphasize that z/WjV 2 ^ is the interaction of vorticity and 
strain since (e.g., Ref. [13J) uV 2 oj= 1/p V X F s , where F? = 2vd / dxk{sik] and p is the 
fluid density. 

In summary, we analyzed small scale enstrophy and strain dynamics in proximity of a 
turbulent /non-turbulent interface without strong mean shear. The experimental results are 
in good agreement with the simulation, at least on a qualitative level, which is considered 
as a clear indication for the reliability of both methods. The behavior of vorticity-related 
quantities is very different from the strain-related counterparts. For example, the viscous 
term is not responsible for building up strain as strain is destroyed by usijV 2 Sij in all three 
regions. In addition, the analysis of these quantities with respect to the distance from the 
interface reveals the range of influence of z/cjjV 2 ^ and z/SjjV 2 Sjj into the non-turbulent 
region. We also found that both UiUjjSij and z/cjjV 2 ^ are responsible for the increase of 
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ij 1 at the interface and substantiate the physical interpretation of the term vujiW 2 uJi as 
viscous interaction, in analogy to UiUJjSij, commonly referred to as the inviscid interaction 
of vorticity and strain. 
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